Translocation of stiff polymers through a nanopore driven by binding particles 
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We investigate the translocation of stiff polymers in the presence of binding particles through a 
nanopore by two-dimensional Langevin dynamics simulations. We find that the mean translocation 
time shows a minimum as a function of the binding energy e and the particle concentration <j>, due to 
the interplay of the force from binding and the frictional force. Particularly, for the strong binding 
the translocation proceeds with a decreasing translocation velocity induced by a significant increase 
of the frictional force. In addition, both e and <j> have an notable impact on the distribution of the 
translocation time. With increasing e and 0, it undergoes a transition from an asymmetric and 
broad distribution under the weak binding to a nearly Gaussian one under the strong binding, and 
its width becomes gradually narrower. 

PACS numbers: 



I. INTRODUCTION 

The translocation of proteins and nucleic acids through 
a nanopore is of essential importance to life, repre- 
sentative examples including the passage of messenger 
RNA through nuclear pores, post-translational trans- 
port of proteins across the endoplasmic reticulum mem- 
brane, and virus injection Jl|. In a seminal experi- 
ment, Kasianowicz et al. [2J have shown that single- 
stranded DNA and RNA can transverse the water-filled 
a-hemolysin channel, which is signaled by a blockade in 
the channel ionic current. Since then, polymer transloca- 
tion through nanopores has garnered high-profile atten- 
tion for its far-reaching technological potential, such as 
rapid DNA sequencing, gene therapy and controlled drug 
delivery 

In addition to its biological and technological rele- 
vance, polymer translocation is also an important issue 
purely from the view of polymer physics. Compared to an 
unconstrained case, the passage of a polymer through a 
nanopore greatly reduces its degrees of freedom and thus 
requires a force to overcome the energy barrier. One 
of the main forms of such driving forces both in vivo 
and in vitro is the trans-membrane force. In biological 
cells, it comes from the trans-membrane electrical po- 
tential, while this is achieved by an electric field mainly 
falling off the pore in experiments, making use of the 
simple fact that biopolymers, such as ssDNA, are nega- 
tively charged. This particular type of translocation has 
been investigated extensively through experiments [U, 
[|[ and theoretical works 0-[I3 ■ Another mechanism uti- 
lizes a chemical potential gradient across the membrane. 
Examples include translocation of chains under different 
solvent conditions 18 20], and that under the chemical 
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potential gradient due to the binding particles (BPs) on 
the two sides of the membrane (9l.[2l|-|42|. In what follows, 
we concentrate on polymer translocation in the presence 
of BPs. Particularly, the chemical potential gradient in 
this work is induced by the binding particles which only 
exists on the trans side of the membrane, as shown in 
Fig.ffl 

As to the translocation mechanism in the presence of 
BPs, two model have been proposed. Simon et al. (2fjj 
suggested that the translocation of chains is a simple 
thermal ratchet process, i.e., the role of BPs bound to 
the translocating chain is only to prevent it from moving 
backward, called the Brownian ratcheting mechanism. 
Later, based on the results from Brownian molecular dy- 
namics simulations, Zandi et al. [3l| considered that the 
binding of BPs onto the chain can provide a mechani- 
cal force capable of pulling the chain through the pore, 
namely, the translocation is a force-driven process. Ad- 
ditionally, some important aspects associated with this 
problem have been investigated successively, including 
the "parking lot effect" due to the size difference be- 
tween BPs and chain monomers [3^ - [39| . the sequence 
dependence (37rl39| , and the chain flexibility [42| . 

Although the studies mentioned above have provided 
a plenty of creative insights into the translocation driven 
by BPs, the underlying translocation dynamics even for 
a stiff polymer still remains unclear. In previous simula- 
tions by Zandi et al. [3l| . the chain length used is fixed 
at N — 16, which is too short, and the binding energy is 
fixed at 5ksT. Therefore, it is very necessary to under- 
stand the influence of the binding energy e, the particle 
concentration cf>, as well as the chain length N on the 
translocation dynamics. 

To this end, we investigate the dynamics of a stiff poly- 
mer translocation through a nanopore in the presence of 
BPs using Langevin dynamics simulations. In section 
II, We briefly elaborate our model and the simulation 
technique. In section III, we present our results and cor- 
responding discussions. Finally, we give a summary in 
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the chain is treated as a completely straight one so that 
we neglect the bending of the chain due to the binding 
of particles. In addition, the pore is considered to be 
completely inert, namely, having no impact on the chain 
other than to permit its one-dimensional fluctuations in 
and out of the cell. 

In the Langevin dynamics simulations, each mobile 
particle is subjected to conservative, frictional and ran- 
dom forces, respectively: 



mfi = -VUi- £v.; 



(3) 



FIG. 1: (Color online) Schematic view of a stiff chain translo- 
cating through a nanopore in presence of binding particles de- 
picted by blue particles inside the cell. The width of the pore 
is w — 1.6a. The radius of the simulative cell is R — 73.5a. 

section IV. 



II. MODEL AND METHODS 

The model geometry we consider throughout this work 
is illustrated in Fig. [TJ where a stiff polymer chain, mod- 
eled as a bead-spring chain of Lennard- Jones (LJ) parti- 
cles with the Finite Extension Nonlinear Elastic (FENE) 
potential, is translocating through a nanopore in the 
presence of BPs. Covalently bonded monomers between 
nearest neighbor beads along the chain interact via the 
FENE potential: 

U FE N E {r) = ~kR 2 ln(l - r 2 /R 2 ), (1) 

where r is the distance between consecutive monomers, 
k is the spring constant and Ro is the maximum allowed 
separation between connected monomers. We employ a 
short-range repulsive LJ potential to incorporate the ex- 
cluded volume interactions between chain segments: 

Here, a is the diameter of a bead, £o is the depth of the 
potential, and r is the distance between monomers. 

Mobile BPs with the same size a as chain segments 
arc modeled as particles freely diffusing within a circu- 
lar cell of the radius R = 73.5a with a pore of width 
w = 1.6a, and they repel each other during the simula- 
tion via the same short-range repulsive LJ potential (Eq. 
2). The circular cell, consisting of one monolayer sta- 
tionary LJ particles, is repulsive for BPs to prevent their 
overflow. The affinity between chain monomers and BPs 
is achieved by exerting an attractive LJ potential with a 
cutoff 2.5cr and the binding energy e. The particle con- 
centration (f> is defined as the particle number density in 
the cell, 4> = N bp /{irR 2 ). Note that, in our simulations, 



Here, U t = Ulj + Ufene(i - l,i,£ + 1) for chain 

monomers and Ui = ^lj f° r BPs, m is the par- 

ticle's mass, which is assumed to be the same for the 
monomer and the binding particle, £ is the particle's 
frictional coefficient, Vj is the particle's velocity, and 
is the random force which satisfies the fluctuation- 
dissipation theorem (43|. The system energy, length and 
mass scales are determined by the LJ parameters So, a 
and bead mass m, leading to the corresponding time scale 
tLj = (ma 2 /eo) 1 / 2 and force scale Eo/a, which are or- 
der of ps and pN, respectively. The reduced parameters 
for all simulations in the present work are chosen to be 
R = 1.5, k = 30, £ = 0.7 and T = 1.2. Thosechosen pa- 
rameters give rise to an effective bond length t of around 
0.96cr. 

The Langevin equation is integrated in time by the 
method proposed by Ermak and Buckholz [44| in one 
dimension for the chain, and in two dimensions for the 
binding particles. This integral scheme was also used 
by Zandi et al. [3lj | . In this way, fluctuations of chain 
monomers along the y axis are forbidden so that the poly- 
mer can only move along the x axis, and such a polymer 
chain is considered to be a rod-like chain. This inte- 
gral scheme means that we project any force on the rod 
along the x-axis and just neglect the y-componcnt. The 
approximation of neglecting y-components of the forces 
does not change the results. Moreover, if the pore is long 
enough, the transverse motion of the rod is not allowed. 

Initially, the first monomer is placed just at the enter 
of the pore (x = 0) and is kept fixed, and the rest of the 
chain and BPs are under thermal collisions described by 
the Langevin thermostat to reach the equilibrium state 
of the system. Then, the first monomer is released and 
both the forward and backward movements of the chain 
are possible. If the first monomer returns to the cis side, 
it would be regarded as an unsuccessful trial. A trial is 
considered to successful as the last monomer exits the 
pore (a; = 0.5 for the first time) and such a trial is called 
translocation. The translocation time r is defined as the 
time interval between the release of the first monomer 
and the exit of the last monomer. Typically, wc perform 
1000 successful runs to get an ensemble average in order 
to reduce statistical errors. 
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FIG. 2: (Color online) Mean first passage time r(s) as a func- 
tion of the translocation coordinate s. The chain length is 
TV = 128. Mean first passage times for the simple diffusion 
and the perfect ratcheting motion are plotted as the dotted 
line (black) and the short dashed line (red), respectively. 
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FIG. 3: (Color online) The mean translocation time as a func- 
tion of the binding energy for three different particle concen- 
trations (j> =2.4%, 12.9%, and 26.8%. The inset is a zoom of 
the data for large e. The chain length here is N = 64. The 
solid lines are guides to the eye, and almost all the errors of 
the data are smaller than the sizes of symbols in the plot. 



III. RESULTS AND DISCUSSIONS 

In our simulations, wc find that the translocation prob- 
ability Ptran, which is defined as the fraction of runs 
leading to successful translocation at given conditions, in- 
creases with increasing the binding energy e and the par- 
ticle concentration </>. For example, at a given (f> = 0.8%, 
we observe P tmn =0.02, 0.20 and 0.44 for e = 1.5, 4.5 and 
7.5, respectively. In addition, at a fixed e = 1.5, Ptran 
increases from 0.02 for = 0.8% to 0.14 for </> = 19.2%. 

A. Mean first passage time r(s) as a function of the 
translocation coordinate s 

We define t(s) as the mean first passage time for s 
monomers exiting the pore. Fig. [5] show r(s) for chain 
length N = 128, binding energy e = 4.5 and different 
particle concentrations <f>. To make a comparison, we also 
calculate the times needed for s monomers exiting the 
pore by the simple diffusion and the perfect ratcheting 
motion, denoted as Tdi//(s) and T ratc het{s), respectively. 
Note that, due to the initialization condition, r(s = 1) 
is the time for the first monomer to first arrive the pore 
exit. Thus, Tdiff{s) = \-°- 5+ ^~ 1 ) e \ (the black dotted 

line), and T ratc het(s) = 0,5 + 2jj~ 1 ^ (the re d dashed line), 
with D = !$jlL being the diffusion coefficient of the whole 
chain. 

As expected, Tdiff(s) is significantly longer than the 
mean first passage time r(s) for different particle con- 
centrations. At a low particle concentration (<f> = 1.5%), 
the translocation is slower than the perfect ratcheting 
motion. However, with increasing c/> from 1.5% to 26.8%, 
the whole curve lies under T ra t c het( s )i demonstrating that 
the translocation is faster than the perfect ratcheting mo- 



tion. This behavior is a little different from the theoreti- 
cal prediction where the perfect ratcheting result is con- 
sidered as the lower limit to the translocation time [35| . 
The reason may be from the non-equilibrium effect of 
the translocation process. We should point out that our 
results are i n ag reement with the numerical findings by 
Zandi et al. [3l[ , where the translocation in the presence 
of BPs is further attributed to a force-driven process. 

Interestingly, at a moderate <f> = 5.4%, the transloca- 
tion is faster than the the perfect ratcheting motion at 
first, but the opposite is the case after s > 85. This dy- 
namical behavior indicates a striking chain length depen- 
dence in this issue. Zandi et al. [3l| have not observed 
this phenomenon possibly due to the too short chain used 
in their simulations (N = 16). 

B. Influence of e and tf> on the translocation time 

Next, we examine the influence of e and <j> on the 
translocation dynamics by measuring the mean translo- 
cation time (r). Fig. [3] shows that for low, moderate and 
high 0, (r) decreases rapidly at first, and then slowly 
approaches a minimum with increasing e. Afterwards, 
a slight but no n- negligible increase is observed, see the 
inset of Fig. [3l Furthermore, the decay rate gradually 
diminishes before reaching the minimum. A similar effect 
of cj> on (r) is observed as shown in Fig. @] 

To understand these behaviors, we first consider the 
relevant timescales and the force from binding. There 
are three relevant timescales: the time To = N8 2 /2D$ 
for a chain to diffuse a distance of a binding site length 
S, with Dq = ksT/t; being the diffusion coefficient of 
a monomer; the typical time r unocc that a binding site 
keeps vacant, and the characteristic time t occ that a 
binding site stays occupied. The typical time T unocc ~ 
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FIG. 4: (Color online) The mean translocation time as a func- 
tion of the particle concentration for three different binding 
energies e =1.5, 4.5, and 7.5. The chain length here is N = 64. 
The solid lines are guides to the eye, and almost all the errors 
of the data are smaller than the sizes of symbols in the plot. 



R b P / D bp = V(0A,p), where R bp 



_1 / 2 is the average 
distance between BPs in solution in two dimensions and 
Db p is the diffusion coefficient of BPs. Taking into ac- 
count the same size for a binding particle and a monomer, 
D(, p . The relationship between T unocc and 

= KT unocc ~ j-^exp(£ /k B T)/D bp [13. 

K = 



we have Dq = 
i~occ meets Tocc 
Here 



%-T-K en is a relevant measure of the cf- 

<j>A c 1 



fective binding strength, K eq = Aoexpfe/kgT) is the 
equilibrium binding constant with Aq = 7rer 2 /4 being 
the typical binding particle area. Then, three dynam- 
ical regimes were divided according to three relevant 
timescales [HI, EE]: the diffusive regime for slow bind- 
ing (to <C t occ , T unocc ), the reversible binding regime for 
fast binding and unbinding (tq 3> t occ , T unocc ), and the 
irreversible binding regime for fast binding but slow un- 
binding (r„„ occ < t <C t occ ). 

For the reversible binding, by calculatin g th e binding 
partition function, Ambjornsson et al. [34l |35| have ob- 
tained the force Fb (in units of ksT/a) exerted on the 
chain by BPs. Considering the univalent binding and the 
coopcrativity effects between BPs bound to the chain, 
for large translocation coordinate s, the finite size cor- 
rections become irrelevant and Fb remains a constant 
value. Fb is given as 



F b {k) «ln{ 



1 



[(" 



) 2 + ( w -l)«] 1 /2}, (4) 



Here, w is a cooperativity parameter. In the present 
work, since BPs bound to the chain interact repulsively, 
0<cj<lifK<f/4, and 1/2 < u> < I for all k values 
[35] . Evidently, k increases with (j) an d £, and thus Eq. 
dU indicates that Fb increases with the binding energy e 
and the particle concentration <J). This is the reason that 
Ptran increases with increasing e and <f>. For irreversible 
binding, however, BPs do not have time to unbind dur- 
ing the translocation, which leads to the thermodynamic 
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FIG. 5: (Color online) Evolution of the polymer flux through 
the nanopore f(t) with the time normalized by the mean value 
for (a) the weak and strong binding, and (b) two different 
binding strengths in the strong binding regime. The inset is 
a zoom of the data in the intersection area of the two curves. 
The solid lines are guides to the eye. The chain length here 
is N = 64. 



evaluation of the force inapplicable. 

Our results shown in Figs. [3] and 0] can be qualitatively 
explained from the perspective of the forces acting on the 
chain during the translocation, including the force from 
binding Fg and the factional force Ff r i c . The initially 
sharp descent of (r) with increasing e and <j> is dominated 
by a substantial increase in Fg. However, the increases 
in e and <f> also simultaneously lead to a larger t occ . As 
a result, Ff r i C rises because segments with bound BPs 
have a higher friction coefficient than free segments (2(| , 
which partially counteracts the favorable factor for the 
translocation caused by the increasing Fg- This fact is 
also the origin of the gradual reduction in the decay rate. 

To further understand why (r) has a minimum as a 
function of e and <f>, we have calculated the polymer flux 
through the nanopore f(t) during the translocation Here, 

f(t) = d ^ t *^ , defined as the differentiation of the average 
translocation coordinate (s(t)) with respect to the time, 
characterizes the instantaneous speed of the transloca- 
tion process, namely, the number of segments passing 
through the pore per unit time. 
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When the binding is weak, it is reasonable to deem that 
the effect of BPs' binding on the friction coefficient of the 
chain is negligible because the unbinding of bound BPs 
is quite fast. Therefore, the translocation under this case 
would be in steady state. However, as the binding gets 
too strong, Ff r i C becomes the prominent factor, leading 
to a gradually decreasing translocation velocity, as seen 
in Fig. [SJa). 

With the two competing effects caused by increasing 
e and 4> operative, we speculate that too strong binding 
is in favor of the translocation in its primary stage, and 
it then results in a more significant increase in Ff ric and 
the ensuing more dramatic decrease in the translocation 
velocity. Fig. (5jb) shows that at a given 4> = 2.4%, the 
polymer flux of the translocation under e = 8.5 is in- 
deed larger than that of the weaker one (e = 6.0) at the 
beginning, but with a faster decay rate. Finally, about 
95% of the translocation under the stronger binding pro- 
ceeds with a slower velocity, and consequently a longer 
translocation time is required. 

Zandi et al. [3l| showed that the mean translocation 
time decreases monotonically with the particle concen- 
tration increasing from 1.1% to 44.2%. In contrast, our 
results indicate that the mean translocation time has a 
minimum as a function of </> which increases from 0.8% 
to 44.2%. The difference is due to the too short chain 
length used in their simulations. As a result, the translo- 
cation has not entered into the frictional force-dominant 
regime. 



C. Distribution of the translocation time 

Previous studies have shown that translocation driven 
by a transmembrane electric field [l?], IH, H|| or a con- 
stant pulling force exerted on the first monomer {47J has a 
qualitatively different shape of the distribution of translo- 
cation time compared to that of the unbiased transloca- 
tion case 0, |48[ . Now that the polymer translocation 
in the presence of BPs is a force-driven process and the 
driving force Fb is a function of the binding energy e 
and the particle concentration </>, it is natural to wonder 
whether changes in e and <f> will have similar effects on 
the distribution of the translocation time. 

As sketched in Fig. H^a), the distribution for e = 2.0 
is asymmetric and broader, while the histogram obeys 
nearly Gaussian distribution as £ increases. Fig. EJb) 
shows that at a relatively large binding energy e = 4.5, 
almost all of the histograms approach nearly Gaussian 
distributions, even under a fairly low <j) = 0.8%. And 
the width of the distribution becomes narrower as the 
binding gets stronger. These results can be ascribed to 
the simple fact that Fb increases with increasing e and 
<t>, and are qualitatively similar to the observations re- 
ported by previous studies [13, IH, H?}- What's more, 
the distribution of the translocation time is more sensi- 
tive to e than <f> in that £ has a more significant effect on 
the effective binding strength n. At a smaller e = 2.0, 




FIG. 6: (Color online) The influence of (a) the binding energy 
e and (b) the particle concentration cj> on the distribution of 
1000 translocation time for a chain of length N = 64 normal- 
ized by their mean value (r). 



we could still observe the transition from an asymmetric 
distribution to a nearly Gaussian one as 4> increases. 

Most recently, it has been demonstrated that the prob- 
ability density function (PDF) of the translocation time 
is solely determined by the dimensionless Peclet num- 
ber, denoted as Pe, [391 ]. It is a dimensionless parame- 
ter comparing drift strength and diffusivity. The PDF 
changes from a broad distribution in the diffusion dom- 
inated regime at small Pe to a significantly narrower 
distribution in the regime of drift domination at large 
Pe. We calculate the effective Peclet number, Pe(<fi c ,e) 
as follows, Pe(2.4%,2.0) ps 8.66, Pe(2.4%,4.5) ps 
15.30, Pe(2.4%,7.5) ps 23.38, Pe(0.8%,4.5) ps 11.55, 
Pe(5.4%,4.5) ps 18.18, and Pe(15.3%, 4.5) ps 21.83. Ob- 
viously, Pe(2.4%, 2.0) ps 8.66 for the broader distribution 
is smaller than these for significantly narrower distribu- 
tions. Therefore, our results about the distribution of the 
translocation time are qualitatively in agreement with 
the findings in Ref. [39|. 
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IV. CONCLUSIONS 

In this work, we have investigated the translocation 
of stiff polymers through a nanoporc in the presence of 
BPs by performing 2D Langcvin dynamics simulations. 
We show that under a certain e and (f>, the translocation 
is faster than the perfect ratcheting motion at first, but 
the opposite is the case for the late stage of the translo- 
cation process. This indicates the striking chain length 
dependence of the translocation dynamics. By scanning 
a large range of the parameter spaces, we find that there 
exist an optimal e and an optimal <j> for the transloca- 
tion. Then, a qualitative explanation is given from the 
perspective of the forces acting on the chain. Finally, 



both e and <p have an notable impact on the distribution 
of the translocation time. Specifically, with increasing e 
and </>, it undergoes a transition from an asymmetric and 
broad distribution under the weak binding to a nearly 
Gaussian one under the strong binding, and its width 
becomes gradually narrower. 
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